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Nomenclature 

Definition 

speed of sound 
wing chord 

section lift coefficient 
wing mean aerodynamic chord 
pressure coefficient 
enthalpy/Q^ 

computational mesh index in the z direction 
maximum value of i-l 

computational mesh index in the r direction 
maximum value of j-1 

computational mesh index in the (f) direction 

maximum value of k in the first quadrant 

index of mesh plane in the first quadrant » neighboring 

the plane ^ y 

index of mesh plane in the third quadrant, neighboring 
the plane 0 « 
llach number 
pressure/p^^ 

perturbation velocity vector/Q^ 

velocity vector/Q 

flow speed/Q^^ 

flow speed 

radial coordlnate/c 

radial coordinate of slipstream boundary/c 
perturbation velocity component in the z direction/Q^ 
velocity component in the z direction/Q^ 
deviation of U from its average value/Q^^ 
perturbation velocity component in the r direct ion 
velocity component in the r dlrection/Q^ 
perturbation velocity component in the 0 dlrection/Q^^ 
velocity component in the 0 direction/Q^ 
axial coordinate/c 
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Nomenclature (Cont.) 


Symbol 

«8 

Y 

5 

n 

n* 

0 

p 

^0 

P 

/s 

P 

4> 

'I' 


Definition 

slipstream swirl angle 
ratio of specific heats 
vorticity vector/ (Q^„/c) 

deviation of slipstream boundary from its 
undisturbed position/ c 
% of semispan 
angular coordinate 

r ^ 1 

“ max 0, (1 “ ^ )J 

stretching factor in r direction 
stretching factor in z direction 
stretching factor in 6 direction 
density/p^ 
density 

perturbation velocity potential/cQ^ 
undisturbed radius of a streaniline/cQ^^ 


Subscripts 

j, denotes inner region variables 

^ denotes out^^r region variables 

s slipstream boundary 

\ denotes total conditions 

oo denotes undisturbed conditions 

Superscripts 

("”) average quantity 






1 . 


Introduction 


Interest In turbo-prop propulsion systems has been recently revived 
due to the propulsion efficiency of these systems and the predicted fuel 
shortages in the future. Propeller slipstreams will Interact with wings 
causing changes in their aerodynamic properties. The deter-inination of 
the effects of slipstream-wing Interaction is therefore necessary before 
any design decisions regarding the installation of turbo-prop propulsion 
systems can be derived. Since current passenger flight cruise speeds 
are in the transonic range, this is the range considered here. 

The interaction of wings with regions of high energy flow embedded 
in uniform streams was studied by Shollenberger (1975), who used flow 
singularities to simulate a Jet interacting with a wing. The flow 
solution and jet position were found by calculating the singularity 
strengths and locations. Lan (1975) used a quasi vortex-lattice method 
and a two-vortex-sheet representation of the slipstream to study the 
interaction problem with different slipstream and freestream Mach numbers. 
Ting, et al., (1972) used the method of asymptotic expansions to study 
the interference of a wing with multipropellers. The effect of nonuniform 
streams on the aerodynamic characteristics of wings nas been studied by 
Chow, etal., (1970) and Kleinstein and Liu (1972) . Boctor, et al., 

(1978) have recent3.y studied the interaction of the slipstream with a 
wing-body configuration at high subsonic Mad; numbers, using the panel 
mt'thod, rnis method does not include treatment of local patches of 

supersonic flow and shock waves. Therefore, though it may indicate 
qualitative trends, it fails to give quantitative descriptions of flows 
in the transonic regime. In addition to the theoretical studies mentioned 
above, experimental tests were conducted by Welge and Crowder (1978) to 
assess the magnitude of the aerodynamic interference of a propeller 
slipstream on a supercritical wing. 

Except for the final experimental study, all previous studies 
dealing with 6lipstream-v;lng Interaction are limited to subsonic flows. 

In this report a simplified model is use.d to describe the Interaction 
between a propeller slipstream and a wing in the transonic regime. The 
undisturbed slipstream boundary is assumed to coincide with an infinite 
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clrculK cylinder. The undisturbed ellpatream velocity le rotational 

turbatlon caused by Introduclns a wine. Into the -Up-t- > also 
rotational. By ™klng small disturbance assumptions, however, the 
perturbation velocity becomes nearly potential, and an ^1”* 
the flow is obtained by solving a potential equation. Th.a s 
::el allows us to obtain basic Information about the Interac . 
problem while avoiding the need tor solving the s»re complex 

equations. 




:.l: 








2. Fonnulation and Governing Equations 

The present analysis considers a wing in a flow with a slipstream 
(inner flow region) embeaded in a freestream (outer flow region). Far 
upstream of the wing, flow conditions are undisturbed by the wing. 

There, the freestream Mach number is , and the slipstream is a 

circular cylinder of radius with velocity distribution 

9iw “ ^ total enthalpy 

velocities and the enthalpy are normalized here by the freestream flow 
speed and the square of the freestream flow speed, respectively, /fhe 
cylindrical coordinates (z, r, 0) are used, with the z--axi8 coinciding 
with the axis of the undisturbed slipstream. The coordinates r and z 
are normalized by the wing chord. The subscripts i, o, <» , and t are 
used 

used to denote inner region properties, outer region properties, 
undisturbed conditions, and total conditions, respectively. As the wing 
is approached, the flow is perturbed from its basic undisturbed condition 
Since the flow is potential in the outer flow region, the velocity 
perturbation there may be expressed as “ V(J>^ . In the slipstream, 
the flow is rotational in general. Let q^ « (u^ , v^ , w^) be the 
velocity perturbation there. Due to the wing effect, the slipstream 
surface deviates from the undisturbed circular cylinder and is defined 
by 

r « R(z,0) « R^ + n(z»0) , 

where T) is the amount of deviation from the undisturbed condition. 

The conditions at the interface specifying zero norraal flow to it 

are 



V. « U.ri 
i i z 
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and the condition specifying zero pressure jump across the interface is 

“ Qi« 2 2 r 

*» H a . • V ' 

2 . o«> i» 


where Q is the magnitude of the velocity Q and a is the sound 
speed normalized by the freestream flow speed. 

The continuity equation governs the flow in the outer region, 
while the continuity equation in addition to the Euler equations governs 
the flow in the inner region in general. The continuity equation is 

V • (PQ) “ 0 , 

or (pW + -H: (■^PV) + i le (P“) - « . 

where p » the density normalized by the freestream density, is given by 

p(z,r,e) - j 1 + ^ [ 1 - ^ ^^2^’f ][ 




ll” 


where 




ai('i') 



* 'i> < K 
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. 'i>> K 

QOO 



. 'i> < Rco 

tl 

, lb > R 


, lb < R 

QOO 

, ilj > R 
♦ 1 00 
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and the unperturbed pressure is given by the relation 


■ i«> 2 

■ ym 

’ o«> 


(which satisfies the Euler equ, tions) , where 

t YW?„(r) 


/ 

R 

CO 


2/ X 

ra^(r) 


dr 


Here, y is the ratio of specific heats, i|;(z, r, 0) is the radius of 
the streamline passing through point ( 2 , r, 0) in the undisturbed flow 
region far upstream of the wing, and the pressure has been normalized by 
the freestream density multiplied by the square of the freestream flow 
speed. 

In general, tbp velocity perturbation inside the slipstream is 
rotational, and a complete solution requires solving the continuity 
equation (Equation 4) in addition to the Euler equation 

(Q • V) Q « - i Vp , r < R , (6) 

Inside the slipstream, solving the continuity equation outside the 
slipstream, and determining the slipstream boundary R(z, 6) . 

2.i Siaall Perturbation Approximation 

Two simplifying assumptions, under which the requirement for solving 
the Euler equation in the inner region is dropped, are introduced. The 
first assumption is that the undisturbed inner region flow is nearly 
uniform: 






^ioo 


where 


U^„(r) - AU^„(r) 






and U„ is the average value of Uj,„(r) . 'fhe second assumption is 

the small disturbance assumption which limits the wings considered here 
to wings with small surface slopes. To see the simplifyir.g effects of 
these assumptions, rake the curl of Equation (6). This gives 

(Q • V)C ® (Q * ^)P * ^^8 ■*" ~2 ^ ^ 

for r < R , where t, is the vorticity and 

^ « V X Q . 

Under the two assumptions made above, Equation (7) to the lowest order 
becomes 


3C 

3z 


0 . 


This allows the velocity perturbation in the inner region to be expressed 
In terms of a velocity potential, and the governing equations simplify 

to ' 


9_ 

8z 


[pW„ + ♦,)] + 7 I? 

7 le -"i v] ■ “ • 


where 


, r < R^ 
“ , r >Il 


W. 


|W (r) , r < 

(O , r > 


and p is given by Equation (5) with replaced by r . The governing 

equation is further simplified by replacing the full potential equation 
with the transonic small disturbance equation 

[p (u, + <!>^) ] + 7 I 7 7 le (r P«>'*’0 ) ” 


dz 


0 , (8) 


where p is to be approximcited by the first four tenns of the blnoiDial 
expansion for Equation (5)* The conditions at the slipstreatn boundary 
reduce tc 


and 


^co - 

^(z, + e, e) 


(R ) 

i<X> ' oo' 


0 , 


<J>^(2, - e. 0) 

R„+ e, 0) “ UiJRoo) 


e 0 


(9) 


( 10 ) 




3. K uJierical Approach - 

The solution to*"the governing equations Is found nuiuericaliy . 

11k 

Therefore, define a net of discrete points (z , r , 6 ) in the c 
putational domain (see Figure 1) with 


0 < i < I 1 , 

1 1 j 1 ^ » 

1 < k < K . 


The computational domain is bounded by an outet cylindrical boundary 
(r “ 9 an upstream vertical plane (z = z ) , and a downstream 

vertical plane (z « . The undisturbed cylindrical slipstream 

boundary (r = R ) lies between the two cylindrical mesh surfaces 
r « r^ ^ and r “ r^ , j “ jg > such that 



'rhe win« extends in both the 8 = ” plane (which bisects the angle between 

It ^ lr4-l 3TT 

the two mesh planes 0 = 6 and 6 ■= 0 , k = k^) and the 9 - 

plane (which oisects the angle etween the. two mesh planes 0 “ 9 and 

e = , k = k^,) . 

3*1 Finite D i fference Formulas 

The finite difference approximation to the governing equation 
(Equation 8) at the point (z^, r^ , 0*^) , 1 £ i 1 I v ^ ^ ’ 

1 < k < K , is 


(gl-t^.j.k ^i44,.j,k^ _ (~l-4.3>k 



where 


» W^(r^) + , 

and denotes an approximation to the function f(z^, , 0 ) . 

The discretized density is first evaluated at points (x — , r , 6 ) . 
The value of is then calculated using the relation 

pi.Jfs.k V ^|pi-!ij,k^pi-»i.^ 


The modified density, , is given by 


gi-t^.j.k ^ ^ ^pi+»S.j,k 


with 


Ap 




(12) 


where 





i-f2 

y 

i+2 

y 


0 

* 0 



and the switching function ]s is given by 


p « max Cl " ^ ) J ' 

The evaluation of the modified density p at points in the constant 

plane, as done here and by vlameson (1976), produces sharoer shocks 
than those produced by the evaluation of p at constant z planes, as 

suGgented by Hafez, et al., (1978). TTie switching function y in 
Equation (12) is evaluated at the point * The evaluation of the 
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swltching function at the point . as done by Jameson (1976). 

produces here a nonstnooth solution at the sonic line. 

Tlie artificial viscosity 


ri.j.k . 


which is an approximation to 






(where Az is the mesh spacing in the z direction), has been added to 
the central difference approximation to Equation (A) in order to produce 
a stable difference scheme in the supersonic zone. This form for the 
viscosity was introduced by Jameson (1976) to solve the full potential 
equation in conservative form. Since the flow is nearly aligned with 
the 2 coordinate, it is sufficient to add the viscosity to the z derivative 

i+^.j,h ’ are 

Tlie perturbation velocity components u and u 

defined by the formulas 


i+1 .1 

z z 


and 




i 1-1 

7.-1 

Sl.iur fomula. for and ara asad ar all »os.. 

points except those ol.cre J - ,i, - 1 »<! J - . respectively. Tl.e 

appropriate for.olas at these points are derived In Appendix A and are 

given by 


i,(J-l)+^,k 


/■vJ 


.J-1^ 


n 4- 


1 ” i 

• J » 


( 13 ) 
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for points in the inner region at the slipstream boundary, and 

.3.3^. 

(r^ - r^“ ) (1 + Ab) 

for points in the outer region at the slipstream boundary, where 


, J - Js » 


A ” u. (R ) 

ioo' a>' 


V.2 2 V 


The perturbation velocity components 
defined by the formulas 


yl.J.kd*! are 


i. j Mh „ JL i_li_ .,.„.uL- — 


W 1 


1 




at all points except those on the wing surface and those in the wing 
wake. On the wing surface the small disturbance boundary condition 


« ^ " ll^(r) / 


is used, where s is the deviation of the wing surface from a horizontal 
pltine and is defined to be positive in the direction of increasing 0 . 
The potential function (f> is discontinuous at the wing wake. There, 
the following formulas are used: 


i.j,k4Ji „ _i , 




k « k 


.3.3A-S . _! . 


k “ k -fl 
w 




1 




. k - k^. 
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^ ^ 4 


6 ^ - 6*^-1 


k ■» k ,-fl 
w* 


where 









» 


and 4^*^ ♦ 4*c ^ are, respectively, the values of the potential ()) at 

U,T) X»,7T 

2 2 

the upper and lower surfaces of the wing trailing edge located at 
0 ~ ^ ♦ r “ r^ , and 7 T respectively the values of the 

’"T '2 

potential 4> at the upper and lower surfaces of the wing trailing edge 
located at <() « , r « r^ • 

The far field boundary conditions are given by 


0,j,k ^ 

u « 0 

and . 1, J+l,k a 

9 ® 0 , 

and on the downstream plane (z « , Equation (11) is solved with 

the z derivative in that equation set equal to zero. 

Finite Difference Fonnulas at the Axis 

Define a Cartesian coordinate system (x, y, z) with a horizontal 
x-axls (0 » y ) and a vertical y-axis (0 » 0) . The continuity 
equation in Cartesian coordinates is 

h [p<“- + »z>] + H <p~'i'x> + h ° • <“> 


The finite difference a|>proximation to Equation (16) is used for points 
along the z-axis and is given by 
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where 


(pv - (pv 


Ax 


(P, 


M )^-»^'^»'^ - (p W 
'^00 c 


00 Q 


Ay 


V » ({) » 

C * 


W « <|) , 

c ^y 


0 , 


( - 1 [( + ( ) 


l.rt.k„.,.l] _ 

( )i>j44s3 „ i j vH-1 ^ j w 


(17) 


and Ax , Ay are the mesh spacings in the x and y directions, respectively. 
The potential function (J) is double valued at the axis* As the axis is 
approached ?'rojn above, It takes a value (t)"^ and as it Is approached 
from below it takes another value ({> • Therefore, set 


,l,l,k 


, 0 < 0j^ < J 

5 ~ < 0 < 2tt 

2 k 


,i- Tr ^ ft c III 
(> » 2 ^ ®k ^ 2^ 


and express the velocity components U , , and in terms of the 

potential function 4> by formulas similar to those used for mesh 


points away from the axis* 
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3.3 Iterative Procedure 

Equation (11) written at all mesh points off the axis and Equation (17) 

written for all mesh points on the axis constitute a nonlinear system of 

ilk 

algebraic equations for the potential function (J> at the mesh 

points. Hafez, et al., (1978) and Holst (1979) have recently solved the 
transonic full potential equation by a simple line relaxation scheme 
which requires no special treatment for supersonic points. This simple 
scheme is used here. Hie finite difference equations at all mesh points 
may be written in the form 

- 0 , J - 1, 2, ... 

, k 1 , 2, « . . 

where is an 1- component vector whose 1^^ element is given by 

the left hand side of Equation (11) if j > 1 and is given by the left 

hand side of Equation (17) if j « 1 . Wlien the velocity components 

appearing explicitly in the finite difference equations are expressed in 

i k 

terms of the velocity potential the vector ^ can be split in the 
following manner 


. J (18) 

, K 



In order to solve Equation (18) iterativelyv a sequence of vectors 

. n « I, 2, ... is defined. Tliese vectors are obti^lned by solving 

I*n 

the linear system of algebraic equations 



, J « 1, 2, . .. 
, k » 1, 2, ... 


. ^1 
. K 
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where is the matrix with its elements evaluated using the 

t* 1 ^ w i k 

iterative solution, and i is the vector with the 

n-l®*" Iterative solution used for evaluating the density and the last 

available Iterative solution used for evaluating the velocity potential. 

The n*^*^ iterative solution is given by 




where w is a relaxation factor. The relaxation sweep is done for 
lines with k”l tok“K on a cylindrical surface and is then 
continued for cylinders with increasing radii. At the end of the relaxation 
sweep, the Jump <n potential at the wing trailing edge 



is found and the n*^'' Iterative value is defined by the relatior 


0^ 

n 


, + 
n- 1 


- "n-l> 


where is a relaxation factor. Relaxation sweeps of the computational 
domain are repeated until convergence occurs. In order to improve 
stability and the rate of convergence, a term 


~ K 


n-1 




<^i- 


Zl_j)At 



1-1, j. 

n-1 



is added in the n*^^ iterative sweep to the left-hand side of the 
algebraic ecjuatlon solved for the mesh point (i,J,h) where K is a 
constant and we set 


At •> z. 


" 1-1 
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A. Results 

The method presented above for solving the slipstream-wing inter- 
action problem has been used to modify a computer code developed by 
M* Hafez (foic solving the transonic full potential equation by the 
artificial compressibility method, described by Hafez, et al#, 1978). 
Numerical examples are presented here for two geometrical configurations. 

The first is that of a simple rectangular wing, while the second is that 
of a swept wing. 

A computational mesh with nommlform spacings in the 6 direction 
and with a stretching factor 

< 6 <1 . 

is used. The mesh is symmetric about both the horizontal and vertical 
planes. Inside a cylindrical region with an axis coinciding with the z- 
axis and with its surface and end planes containing the outermost wing 
tips, the mesh spacings in both the axial and radial directions are 
uniform. Outside this cylinoer, however, the mesh spacings are non- 
uniform. They are stretched with a stretching factor 5 in the 
positive radial direction. Downstream of the wing the mesh spacings 
are stretched by the factor in the positive z direction while they 

are stretched by the same factor upstream of the wing in the negative 
z direction, llierefore the mesh spacings Increase in size as one moves 
away from the wing. A value of 1.1 for each of three stretching factors 
is used here. 

The results are calculated for a freestream Mach number M “0.8 

QOO * * 

The slipstream swirl angle and total pressure distributions 

of highly loaded propellers under development (see Welge and Crowder, 

1978) are shown in Figures 2 and 3. These distributions are used in the 
calculations presented here. It is noted that Equations (13) and (14) 
are valid for general A and B values. However, for the p /p 
distribution of Figure 3, both A and B are approximately unity. In 
this case, it is consistent with the small disturbance approximation to 
set A - B « 1 , allowing the problem to be solved with no special 
treatment for points at the slipstream boundary. Although the application 
of Equations (13) and (14) is a simple matter in the present calculations 
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„her^ cylindrical coordinates are used, a more complex set of equations 
replaces them when general coordinate systems are used. In such cases, 
a great simplification is caused by avoiding any special treatment for 

points at the slipstream boundary. 

The first set of results are calculated for a rectangular «lng »lth 
»ACA 0012 sections, a 3° angle of attack, and a .arlnun swirl angle 
ew » 3° . The rectangular wing planfor.m is sketched in Figure A. 

■rJe^alculatlons for the rectangular wing are carried out on a 49x20x24 
(a t, 0) mesh, with 14 nesh points along the wing chord and 20 po nta 
alLg’the wing span. The conputatlonal mesh extends 1.5 chord lengths 
upsrrea. of the wing leading edge, 5 chord lengths down.treaw of the 
wing trailing edge, end the outer cylindrical computation mash boundary 
radius is 4 chord lengths. The first iterative guess Is taken to e . 
linenrly Interpolated solution obtained by solving the problem on . 

rough mesh 40x16x16 

Profiles of the pressure coefficient 

Cp « 2(p - P^) 

at aeetlons A and B of the wing (see figure 4) are given In Figures 5 
and 6 . At Section A, the swirl angle effectively reduces the wing s 
angle of nttack. On the upper surface, the axl.l velocity lnc«ment lu 

the slipstraam and the swirl angle produce opposite effects. ey, 

however, produce similar effects on the lower surface. The swirl angle 
rends to produce en upstream displacement In the upper surface » “ ^ 

poaltlou. While the axial velocity increment rends to produce a downsttea 
dlaplacemenr. At Section B. the swirl angle sffectivel, increases the 
wing's angle of attack. The axial velocity Increment and the awlrl 
angle produce aimllar effects on the wing's upper surface and opposite 
affects on its lover surface. Note in Figure 6 that f,.e slipstream 
effect produces a downstream displacement in the shock position, and 
pcoduc. a .trongor shock. This cn„se,uontly contributes to eu lucrease 
in tbo wave drag. Figores 7 and 8 show the dlstrlbuti.tn of 
(the mini™ value of C along a chord) along the wlog apsn. In tne.e 
figures, we see that the slipstream effect is largest inside the si p- 
atren. but extends outside It. Figure 9 show, the distribution along 
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tK. vl„s span »£ <=. • ‘"a ''atinn U£t co.£flcla„£. It in t-lnded 

f.„„ this flguta thL th. supstrean effect on lift distribution Is 

primarily due to the swirl angle. . u ^r. 

The effect of the slipstream on the wing's wave drag s 8 own 

leble 1. Estlnateo of ..ve dreg are found (see Muruan and Cole, 
by integrating 

.3 


/ 

shock 


4 >, 


ds 


ubere the Integral la taben along sbocb u.ves and 1*^1 
d, across the shock. In the table, the drag Is nor^lited by that du 
rig seulspau In the absence of the sllpstreau - -- -- 
that the sulrl effect on uavo drag Is that of rodlstclbutl n. 
effect o„ the total »ave drag, hovever. 1, suall. This la In con rast^ 
to the sllpstreau arlal velocity uhleh produces a large Increase In »a 

''“'■■The effect of the sllp.trea. on the ulng's lift Is show. 1„ Table 2. 

Esclles of lift are found by integrating tulce Mee <fa 

In 4. betucen the upper and lover ulng trailing edge Itbat 

the ting span. The lift values show. In Table 2 are nor.all^ at 

due to the Ulng seulspan In the absence of the slipstream. The table 
Chous that the effect of the sllpatreau on the ulng s total 

second example, attempting to calculate a case uhleh Is as close 
as possible Co the erperii.ental configuration of Uelge and Crou er, 

, 1. nou presented. The pla,.fotm of the ulng-body conflgnrat In 

t; Velge aid Ctouder aPlS, la depleted In figure ,0 ue o 
simplicity of the computer code In the present study, It s not posslhl 

to slmnlate flous about complek geometries. The calculat, one aye 

therefore been done for a simple swept Isolated wing conflgnrat on. 

Ulng planform 1. ske-tched in Figure U. and the defining 
airfoil section is given in Table 3. This airfoil section s e . 
as that of the .aperlmentol ulng section at the 35 percent soml.pan 
atatll The nse of a Cartesian coordinate system at the sXfpstrea. 
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axis, in the manner discussed in Section 3.2, has been found to cause 
solution inaccuracies near the leading edge of a swept wing. To avoid 
this problem, an infinite cylinder whose axis coincides with the slip- 
stream axis is assumed. The solid wall boundary condition is applied at 
the cylinder's surface. The cylinder's radius is chosen to be half a 
mesh spacing in the radial direction. A wing angle of attack is 

chosen to roughly match the sectional lift coefficient in the slipstream 
region to that of the experimental lift coefficient in the absence of 
slipstream effects (swirl and axial velocii./ Increment) . In the calcula- 
tions, the total pressure distribution in the slipstream is assumed to 
be that depicted in Figure 3. Calculations with both positive swirl 

angle (a ”7°) and negative swirl angle max “ ^ 

s m£tx 

been performed. 

The calculation for the swept wing are carried out on a 107x20x20 
(z, r, 6) mesh, with 16 mesh points along the wing chord, and 20 points 
along the wing span. The computational mesh extends 2.5 chord lengths 
upstream of the wing leading tip, and 8 chord lengths downstream of the 
wing trailing tip, and the outer cylindrical computational mesh boundary 
radius is 5 chord lengths. The first iterative guess is taken to be a 
linearly Interpolated solution obtained by solving the problem on a 
rough mesh 78x14x14 (z, r, 6) . 

Profiles of the pressure coefficient at sections A, B, C, A’, B', 

and C' of the wing (see Figure 11) are given in Figures 12 through 17. 

A comparison of the solutions at sections A and A’, which are closest to 

the axis, indicates a strong effect due to the cylinder co-centered with 

the axis. The eifects of this cylinder is to compress the flow at 

section A, and to expand the flow at section A'. These effects may be 

concluded from basic flow properties along swept wings and are confirmed 

by comparing the solutions presented in Figures 12 and 15. At sections 

A, B, and C, the positive (negative) swirl effectively produces a reduction 

(increase) in the wing's angle of attack, while the positive (negative) 

swirl effectively produces an increase (reduction) in the wing's angle 

of attack at sections A', B', and C'. Figures 18 and 19 show the distii- 

bution of C , along the wing span, and Figure 20 shows the distribution 
p roin 
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of the sectional lift coefficient along the wing span. The effect of 
the slipstream on the wing’s lift is shown in Table 4. Drag calculations 
have not been performed, since the solutions indicate supersoaic-supersonic 
shocks. The simple calculations which are used to estimate wave drag 
due to supersonic-subsonic shocks arc not applicable in the case of 
supersonic-supersonic shocks. 

Results presented in this example show similar trends (for slip- 
stream effects on supercritical wings) to those of the experimental 
results produced by Welge and Crowder (1978) (see Figures 21 and 22), A 
qualitative comparison between the experimental results and the numerical 
results is not appropriate due to differences in the geometrical configura- 
tions of the experimental and numerical test cases. The absence of a 
body next to the wing and the presence of a cylinder in the slipstream 
in the calculated example, in addition to differences in the wing planforra 
and its cross-sectional profiles between the experimental and numerical 
examples, contribute to differences in the flows about the experimental 
and numerical configurations. In addition to the current code’s geometri- 
cal limitation, it should be noted that the use of the small disturbance 
equations near Lae blunt leading edge of a supercritical airfoil causes 
solution inaccuracies in the leading edge region. Precise numerical 
calculations for complex geometrical configurations may be made by the 
use of computer codes capable of handling these geometries. Leading- 
edge inaccuracies may be avoided by using full potential codes. 

Tlie calculations were done on a CDC 7600 computer. For the rectangu- 
lar wing, 1 second and 2 seconds were required per iteration, respectively, 
on the rough mesh and the final mesh. The number of iterations required 
to reduce the maximum residual to 0.001 was approximately 150 for the 
rough mesh and 100 for the final mesh. For the swept wing, 1,36 seconds 
and 4.6 seconds were required per iteration, respectively, on the rough 
mesh and the final mesh. The number of iterations required to reduce 
the maximum residual to 0.001 was apprc. ^mately 7'; for t!ie rough mesh 
and 190 for the final mesh. Tlie values used for the relaxation factor 
03 varied between 1.4 and 1.6. Tlie relaxation factor Q for calcr.lating 
the jump in potential at the trailing edge was taken to be 1,2, rv;d 
K , the <j)^^ term coefficient, was set equal to 1.0. 
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5. Concluding Remarks 

In this report, an inviscid model for the interaction between a 
thin wing and a nearly uniform propeller slipstream has been presented. 

In this model, the perturbation velocities due to the interaction are 
potential, even though the undisturbed slipstream /eiocity is rotational. 

This allows basic information about the interaction problem to be obtained 

while avoiding the need for solving the Euler equations. Tor typical 

slipstream velocity distributions, only minor modifications to the free 
flow potential equation and wing boundary conditions are required to 
produce the slipstream effect. These modifications can be easily incor- 
porated into available transonic potential codes. 

The slipstream effect on a wing in the transonic regime has been 
demonstrated through a simple example. Solutions obtained for a rectangu- 
lar wing indicate that the slipstream has a strong effect on the aero- 
dynamic properties of the wing portion immersed in the slipstream. The 
effect of the slipstream on the rest of the wing is less, but t is 

discemable. The results indicate that the slipstream swirl has a 

strong effect on the wing load distribution, however, its effect on the 
total wing lift and wave drag is small. 'The axial velocity Increment 
inside the slipstream has little effect on the wing's lift; however, it 
causes a large increase in wave drag. 

Although the computer code used for the present study is a simple 
code which does not alxow calculations for complex geometrical configurations, 
an attempt has been made to calculate a case which is as close as possible 
to the experimental configuration of Welge and Crowder (1978). The 
calculated results indicate the general trends for slipstream effects on 
a supercritical swept wing. However, solutions calculated are not 
expected to be accurate near the blunt leading edge of the supercritical 
wing, since the small disturbance equation is used in the present calculations, 
Therefore, the change in leading edge suction associated with the slipstream 
swirl is not properly accounted for in the solution. The computer code 
used in this study is a simple code for the purpose of demonstrating the 
interaction effects for simple geometrical configurations. More sophisticated 
codes are available (see Jameson and Caughey (1977) and Caughey and 
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Table 3: Airfoil S 




0.0 

0.05 

0.25 

0.50 

1.25 

2.50 

5.00 

7.50 
10.00 

15.00 

20.00 

25.00 

30.00 

35.00 

40.00 

45.00 

50.00 

55.00 

58.50 

65.00 

70.00 

75.00 

80.00 

85.00 

90.00 

92.50 

95.00 

97.50 
100.00 


Upper Surface 
Coordinate, y 


-0.0568 
-0.0540 
-0.0491 
-0.0455 
-0.0393 
-0.0327 
-0.0245 
-0.0189 
-0.0144 
-0.0079 
-0.0033 
0.0000 
0.0023 
0.0035 
0.0041 
0.0035 
0.0022 
0.0001 
- 0.0020 
-0.0069 
- 0.0121 
-0.0632 
-0.0257 
- 0.0344 
-0.0444 
-0.0502 
-0.0565 
-0.0637 
-0.0713 


Lower Surface 
Coordinate, y 


■ 


-0.0568 

-0.0619 

-0.0670 

-0.0708 

-0.0780 

-0.0860 

-0.1256 

-0.1030 

-0.1081 

-0.1151 

-0.1197 

-0.1228 

-0.1250 

-0.1260 

-0.1260 

-0.1248 

-0.1223 

-0.1179 

-0.1138 

-0.1041 

-0.0952 

-0.0860 

-0.0771 

-0.0696 

-0.0646 

-0.0640 

-0.0651 

-0.0682 

-0.0734 

















\ 


- 29 - 


Figures 
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Rgure 1. Cross Section Through Computational Mesh at a Constant z Plane 
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Section B 



Section A 
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Figure 4. Rectangular 



ir^g Planform, 
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Figure 8. Minimum C„ Distribution Along Rectangular Wing Bottom. 









Section C' 
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Uniform Flow 
olipstream with Positive 



Figure 21. Wing Surface Pressure Variation Due to Power and Swirl at Mqoo -0.8 
(Experimental Results Taken From Wedge and Crowder, 1378). 



Figure 22. Span Loading at Mqco =0.8 (Experimental Results Taken From 
Wedge and Crowder, 1978). 
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Appendix A 

In general, the approximation of the radial velocity component at 
the point (z^, 6*^) is given by 

rJ - 

This, however, is not true at mesh points fox which 1 « 1 since the 

s 

slipstream boundary (across which the velocity components are discon— 

tinuous) coincides with the cylindrical surface r * r^ ^ , where j j 

Expressions for the radial velocity components v-^*^ in the inner 

region at the slipstream boundary and in the outer region at the 

slipstream boundary may be found in terms of the potentials and 

i 1~1 k 

<{> * * > where j « , as follows: 

® J 1 

Define a fictitious potential at the point (z^. r^. e*^) , 

where j = , so that 





It follows from Equation (9) that 



(A.l) 



A 3> 


i,k 


j-1 

r J 



(A. 2) 


Tl>e axial velocity perturbation components , in the inner 

teglon at the slipstream boundary, and , in the outer region at 

the slipstream boundary, may now be expressed as follows: 

/(r^ - r^-l) 


rJ - 
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and 


. [♦‘-1.3.'^ . ^ v+l-l-V - rl-h] { 

/ (r^ - r^”^) 

,j J-1 2 ^ 2 ^ 


r- - r 




+ “ 4 > 


r 


) . j - jg . 


Substituting the expressions tor and 11 +-^'’'^’*' into Equation (10) , 

the following relation is obtained: 

+ AB) - + (1 - 

» .(. AB) - 2B()>^~^’^’’^ + (1 - 


. j “ Js • 


(A.3) 


Ni.k 


Tills all ws |J> ’ to be expressed as follows: 


vi.k 


(G + 2b4>^’^’’" - (1 - ab) 4,^*^-^’S . j o . 


Substituting into Equations (A, 1) and (A.2), the radial velocity components 
become 


i,k 


V- 


^i*k 
v4- ’ w 


(rJ - r^“^) (l +AB) 


(r-^ - r^“^) (1 + AB) 


- (G + 2B'|)^*J’‘‘ -2(}>^’^"^’’") , J » j. 
(G + 2B(|)^’^’*^ - 2(|.^’^”^’*^) , j “ j 
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Since relation (A.3) is true for all values of i , it may be 
applied repeatedly for decreasing values of i » obtaining 


G » + AB) - + (1 - 


, J - J, 


Since the radial velocity component is zero at the upstream boundary* 
this leads to the conclusion that 


KO.J.k 


, j “ jg . Jg + 1 . . . j + 1 


Since 




^o,J+l,k ^ ^ 


it is concluded that 


C » « 0 > J “ J 


G « 2c^ ; 


The constant * and therefore G , may be chosen arbitrarily. Choose 

: ■ G « 0 . 

The final expi*essions for the radial velocity components at the slipstream 
boundaries are therefore 



f j “ j„ . 

(i-J _ r-^-1) (i +AB) 

s 


. J “ • 

(r^ - r^"^) (1 + AB) 

s 


./ 
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